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ABSTRACT 

In this paper an automatic method applicable to the segmentation of mediastinum Computerized Axial 
Tomography (CAT) images with tumors, by means of Alternating Sequential Filters (ASFs) of Mathematical Morphology 
and connected components extraction based on continuous topology concepts are described. Digital images can be related 
to topological space structures and then general topology principles can be straightforwardly implemented. This method 
allows not only to accurately determine the area and external boundary of the segmented structures but also to obtain their 
precise location. 
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1. INTRODUCTION 

Technological development has significantly improved diagnostic imaging, enabling renal tumor and incidental 
hepatic tumor detection (usually small in size) in younger people and with an eventually lower malignant potential. This 
has led to a remarkable advance in interventionist techniques such as cryosurgery and radiofrequency ablation, preventing. 
Notwithstanding this, both cryosurgery and radiofrequency ablation, through extremely low and high temperatures 
respectively kill tumor as well as healthy cells, rendering crucial the identification of tumors with an extraordinary spatial 
accuracy. 

The medical image played a pivotal role in the early detection, diagnosis and treatment of cancer. Under certain 
circumstances, early detection through medical imaging allows tumor cure or elimination. The new interventionist 
techniques such as cryosurgery and radiofrequency ablation, in some cases, prevent major surgeries. However, both 
techniques, through extremely low and high temperatures respectively, kill tumor as well as healthy cells, rendering crucial 
for the identification of tumors with an extraordinary spatial accuracy. 

In CAT imaging, different tissues offer no uniformity and contain intensity variations. Conventional image 
segmentation techniques rely on pixels properties treated on an individual basis and hence fail to segment [1],[2],[5]. 

In the first stage, opening and closing Alternating Sequential Filters (ASFs) are applied based on Mathematical 
Morphology Reconstruction [8]. The advantage of this kind of sequential morphologic filters is to filter the desired objects 
leaving the original shape unaltered. This results in less image distortion [12], [13], [15]. ASFs by Reconstruction consists in 
opening and closing iterations by Reconstruction with structuring elements of growing size, being necessary to use a 
second structuring element in the Reconstruction operation. In this way, the regions of the image describing meaningful 
details remain connected [9] . 
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In the second stage, an algorithm is applied along the lines of the connected components criterion on the basis of 
continuous topology. Once the connected components are obtained, the boundary and interior of each of them are 
calculated, so yielding the desired image segmentation. 

This work is organized as follows: Section 2 covers the theoretical principles applied, Section 3 deals with the 
methods suggested, Section 4 presents the results yielded by the different images employed and Section 5 lays down the 
conclusions drawn. 

2. PRELIMINARIES 

This section presents some general definitions and general topological results, which will be used along the 
development of this article. Among the topologies which can be regarded in a set X, there are two representing opposed 
ends. One is the discrete topology J =P(X), in which all parts of X are regarded as open sets, and the other is the indiscrete 
or chaotic topology, J '= {0, X}, in which 0 and X are the only open sets. Indeed, J and J ' represent, respectively, the 
major and minor topologies that can be defined in a set. 

A set U of a topological space ( X , J) is a neighborhood (J-neighborhood) of a point v if and only if U contains an 
open set to which v belongs. U{x) shall denote the neighborhood system of point v. 

2.1.1- Topology Based on Neighborhood Systems 

Let (X, J) be a topological space, for each v eX, let U(x) be the neighborhood system of point v. Then: 

(i) If £/ e U(x) then xcU 

(ii) If U and V are members of U(x ), then U n V c U(x). 

(iii) If U e U(x) and f/cf then V c U(x). 

(iv) If U cU(x), then there is a neighborhood V eV(v), such that Va U and V e U(y) for each y of V (i.e., V 
is neighborhood of each of its points). 

If N is a function assigning a non-empty collection U(x) to each r e X , of parts of X satisfying (i), (ii), (iii), then 
the collection J of those sets U such that U e U(x) for each v of U is a topology of X. If (iv) is also met, then U(x) is 
precisely the neighborhood system of v with respect to topology J[3]. 

2.1.2- Topology based on Distance Function 

There are many topological spaces in which topology are deducted from a distance notion. Distance, defined in a 
set X, is a function d:X x X — ► R such that for all x, y, z eX, the following is verified: 

(a) d(x, y) > 0 

(b) d(x , y) =0 =>v =y (separation) 

(c) d(x , y) =d(y , x) (symmetry) 

(d) d(x , y) <d(x , z) +d(z , y) ( triangle inequality) 

A metrical space is a pair (X, J), such that J is a metrics for X [7]. 


Impact Factor (JCC): 2.6305 


NAAS Rating 3.19 



Topological Digital Image Processing 


55 


Let’s see now how, from a given distance, topology can be generated. If a distance d is defined, it is possible to 
define the following: 

B(x,r ) = { y eX; d(x,y)<r} open ball of center x and radius r. 

By means of these, opens set can be generated: a set A czX is open, if and only if, for every xc A there is an open 
ball of center v and radius r included in A. Then the collection made up of all open sets constitutes a topology termed 
topology associated to distance d [ 7]. Given the fact that the infinite collection of open sets generated is different for each 
distance, consideration should be given to the likelihood of associating a different distance topology to each of them. 

If (X, J) constitutes a topological space, and Y is part of X, it is possible to build a topology J y of Y called relative 
topology or relativization of J to Y. Relative topology Jyis defined as the collection of all member intersections of J with Y. 
Every member U of relative topology J y is regarded as open in Y. The topological space (T, J y ) is called subspace of (X, J) 
space. 

A partition of a set X is defined as its own division in classes X h X 2 , X 3y ... X n such that if x cX , x belongs to one 
and only one of these classes. 


Formally, a partition of X is defined as a function p: X — > J, such that \/x eX, x cp(x) and Vx, y eX 
p(x) =p( y) o p{x) fl p( y) = 0being p{x) the partition of X containing x(See figure 1). 



x 3 



A set SaX is said to be connected if it cannot be written as the union 
of two disjoint non-empty open sets. Another possible way of 



Xi __ 


expressing set connection is by stating that it is possible to go from 
one set point to another one by means of a continuous movement, 


*4 \ " 

x 5 


not leaving the set. This leads to the definition of a path connected 
space, or simply, connected space. 

A path in a topological space X is a continuous function/- [0,1] — >X. 

Figure 1: Example of Partition of a Set X 

Points a =/( 0) and b =/(l) are the path ends/. In this case, the path 
is said to link point a with point b. When a = b, f is said to be a 
closed 


..AO) 

1 ,/(0)=g(0) 

g(l) 


path. 

Given the paths / g : [0,1] — >X with/(l) =g(0) , path fvg : [0,1] — >X 
can be defined as resulting from following path / first, and 
afterwards path g (See figure 2). 


/vg:[0,lHI 



Formally, a set SczX is regarded as path-connected when any two 


Figure 2: Example of Path- Connected 
Spaces 


points of S can be linked by path /in S. 


A topological space X is regarded locally path-connected when for every v eX and every neighborhood U of x 
there exists a path-connected neighborhood V such that x c VdU . 
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Interior, set boundary and isolated point are basic concepts in a topological space. A point x eDcA of a 
topological space is an interior point of D if and only if D is a neighborhood of x, and the set of all interior points of D is 
the interior of D , denoted by: int(Z)) . The interior of a set D can also be defined as the union of all open sets contained in 
D, i.e., int(D) is the major open set in D. 

Set D closure is defined as the minor closed set CcA, such that DczC. The boundary of set D is set: 
d(D) =closure(D) -int(D) 

A point x of D cX is isolated if there is a neighborhood U c U{x ), such that U \ {x} HD =0. 

Given an arbitrary point x cD czX , there are three possibilities mutually excluding: there can either be a 
neighborhood U cU(x) contained in A (i.e., x e int(D) ) , a neighborhood U cU(x) contained in (X -D) , (i.e., x c int(X -D) ) 
or every x neighborhood contains points of D and of (X -D) , (i.e., x c d(D) ). Then, every D set determines space 
decomposition in three disjoints subsets two to two [14]: 

X = closure (D) U int(Z)) U int(X -D) 

2.2.2- Topology in a Discrete Space Z 2 

A digital image is a function f:Z x Z— ►[(),... N - 1] in which A -1 is a positive whole belonging to the natural 
interval [1,256]. Considering Z x Z with topologies defined by metrics cT and d” restricted to Z x Z defined by: 

2 

d'ix.y) \x t -y t \ 

i= 1 

d"(x,y) = max \x t -y t \ 

the following neighborhoods of x EZ x Z can be defined, respectively (See figure 3): 

u 4 (x) = { y e Z 2 / lyl - xll + Iy2 - x2l < 1 } 

= {yeZ 2 /d’(x,y)<l} 

Ug (x) = { y e Z 2 / max [lyl - xll + Iy2 - x2l] < 1 } 

= { y £ Z 2 / d”(x,y) < 1 } 

Generally, neighborhoods characterization is geometrically understood as an “approximation”. 

However, when images in levels of gray are employed, 
neighborhoods depend on them. Neighborhoods can be 
defined by means of intensity, as follows: 

/: G cZ 2 — 2V -1] 

U f (v) =f’(f(v) - c,f(v) + c)nU(v), c e [0,... N-l] 

In a discrete plane, a path is made up of a sequence of 
points 
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Figure 3: Different Neighborhoods of a 
Point x Depending on the Metrics 
Employed 
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a = x 0 ,...,x k = b in which x t lies adjacent to x i+1 fov i = 1,..., k (a y point is adjacent to x, if there is a 
U x neighborhood, such that y eU x ). 

Let XaZ x Z be an image, two points x and y are said to be connected in X if there exists a path included in 
X , joining x with y. If a topological space X is not connected, we can wonder how many “parts” X has. Space parts are its 
connected components. Let X be a topological space and x a point of X, a connected component of x eX is the union C x of 
all connected subsets of X containing an x. 

There is at least a connected subset of X with an x: {x}. Then C x is not empty. C x is the major connected subset of 
X containing an x. 

2.2.3. Grouping Criterion 

Let X be an image in levels of gray and Ja topology associated to X ; and let SO be defined as cp: S x X — >R, such 

that: 

V(A,x)=d(x(A),x) (1) 

Where y(A) is a characterization of A and d a metrics. 

Then cp is said to be a function yielding the set similarity relationship A cS with x eX, after considering y(A). This 
function depends on the application under consideration. 

Given a fixed £ and 5 and SCJ. Let X be an image in levels of gray and AcS, an element xeX is said to belong to A 
if the following is met: 

B(x, £) HA \ {x} +0 y cp (A \ {*},*)< 6 

i.e., S is a covering of X. In other words, for each x cX, there exist AcS such that xcA. 

Topologically, this definition determines that an element in the image will belong to an element AcS if there is a 
ball with a center in x and a given radius, such that interception occurs; and that said element is related to it. 

3. MATERIALS AND METHODS 

In the first stage, reconstruction opening and closing Alternating Sequential Filters (ASFs) were applied by 
reconstruction. AFSs consist in the iteration of morphologic closing and opening operations by reconstruction, with 
structuring elements of growing size, being necessary to use a second structuring element in the Reconstruction operation. 
The advantage of this new type of filters is that they filter the desired objects leaving the original shape unaltered, thus 
resulting in less image distortion. In this way, the regions of the image describing meaningful details remain connected. 

In the second stage, the proposed algorithm employs the grouping criterion defined in the prior section to obtain 
the desired image segmentation. The algorithm starts by generating a set with an image point (the first pixel found at the 
beginning of the image, or that pixel verifying a certain condition) assuming certain order or direction. Based on this, it 
creates the image connected components according to the grouping criterion above defined in equation (1) for a specific 
function cp defined in this article. 

Algorithm to extract connected components of the image: 
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To begin with, function cp should be defined as it determines the grouping criterion, i.e., distance d and y(A) • Let 
d G be the geodesic distance [4] [6] [12] and y (A) = {a c c A: a c is center of A}. 

Given a fixed £, a first set A n = {x} is created, such that x eX, with n-1. 

Step 1: 

If y EX such that y £ Af/I = 1 . . .n, go to steps 2 and 3, 
if there is no yEX go to step 4. 

Step 2: 

For I = 7 . . . n if B(y,£ ) C\A t f 0 and cp (A h y ) <oo, in I -l U { i} (being I a set of indexes). 

Step 3: 

If I f 0 then Aj = AiU{y} V/ El; otherwise 
n = n +1; A n = A n U{y} 

Step 4: 

Once S = Ui EI A j covering is obtained, the connected components Cj are generated. 

Step 5: 

Once the resulting connected components Cj are obtained, the boundary and interior of each of them is obtained. 

Step 6: 

Calculation of the mass center for connected components j-esima. (whose localization information is of interest 
when it comes to practical ends). 

Step 7: 

Visualization. 

The algorithms were programmed using Matlab5.3. Standard functions of this language and a specific library 
called SDC Morphology Toolbox (SDC, 2001) is employed. 

4. RESULTS AND DISCUSSIONS 

Figure 4 lists the results obtained with the proposed method in a CAT image of the mediastinum, in which a liver 
tumor can be observed. First the original image is presented, then image b) after applying the Alternating Sequential Filters 
with structuring elements of growing size. This allows homogenizing the areas to be detected without subjecting them to 
deformations. Image c) depicts different connected components marked in different colors. Finally, image d) shows the 
correct tumor segmentation. The exact boundaries of the segmented tumor as well as its localization can be calculated in 
relation to it. 
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a) Original Image 


b) Image Obtained after 
Applying Sequential Filters 


c) Image with Different 
Connected Components 


d) Image Resulting from 
Liver Tumor Segmentation 


Figure 4: CAT Mediastinum Image 


Figure 5 depicts several CAT mediastinum images and their corresponding segmented images with the proposed 


method. Original images can be observed in the first column; while the second column shows the images resulting from the 


segmentation of the different tumors (lung, kidney and liver, respectively). The third column depicts the contour extraction 


of said tumors. As seen in the figure, the proposed segmentation is genuinely satisfactory in all cases, having been tested in 


a group of more than 50 testing images. 




Images - Original 


Images - Segmentation of Different Tumors 
(Lung, Kidney and Liver, Respectively) 


Images - Perimeter Extraction of the 
Tumors Previously Segmented Can be 
Observed 


Figure 5: Different CAT Images 


Column - 1 


Column - 2 


Column - 3 


5. CONCLUSIONS 

This article describes an automatic method applicable to the segmentation of mediastinum Computerized Axial 
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Tomography (CAT) images with tumors. It has proved to be accurate and efficient, thus enabling final users to know the 
tumor area and localization with remarkable spatial accuracy. This segmentation method turns out to be optimum in tests 
made prior to cryosurgery or radio frequency ablation. 

The potential of this method lies in the possibility of employing different distances and characterizations of 
mounting interest in the segmentation of other types of images. For instance when cp is defined by the Euclidean distance 
and y(A) by a level of gray, the grouping algorithm is that known as image labeling. 
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